Continuous Paranematic-to-Nematic Ordering Transitions of 
Liquid Crystals in Tubular Silica Nanochannels 
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The optical birefringence of rod-like nematogens (7CB, 8CB), imbibed in parallel silica channels 
with 10 nm diameter and 300 micrometer length, is measured and compared to the thermotropic 
bulk behavior. The orientational order of the confined liquid crystals, quantified by the uniaxial 
nematic ordering parameter, evolves continuously between paranematic and nematic states, in con- 
trast to the discontinuous isotropic-to-nematic bulk phase transitions. A Landau-de Gennes model 
reveals that the strength of the orientational ordering fields, imposed by the silica walls, is beyond 
a critical threshold, that separates discontinuous from continuous paranematic-to-nematic behav- 
ior. Quenched disorder effects, attributable to wall irregularities, leave the transition temperatures 
affected only marginally, despite the strong ordering fields in the channels. 

PACS numbers: 61.30.Gd, 42.25.Lc, 64.70.Nd 



Spatial confinement on the micro- and nano-scale can 
affect the physics of liquid crystals (LCs) markedly. Mo- 
dified phase transition behavior has been found in ex- 
jeriments on LCs imbibed into a variety of porous media 
1 , 0] , in aerogels 0] , in semi-confined thin film geometries 
4 1 , and at the free surface of bulk LCs ^ . 

For example, the heat capacity anomaly typical of the 
second-order nematic-to-smectic-A [N-SmA) transition 
in rod-like LCs immersed in aerogels is absent or greatly 
broadened. This allowed a detailed study of the influ- 
ence of quenched disorder introduced by random spatial 
confinement on this archetypical phase transition H. 

It has also been demonstrated experimentally [l|, Q , 
in agreement with expectations from theory 0,01, that 
there is no "true" isotropic- nematic (I-N) transition for 
LCs confined in geometries spatially restricted in at least 
one direction to a few nanometers. The anchoring at the 
confining walls, quantified by a surface field, imposes a 
partial orientational, that is a partially nematic ordering 
of the confined LCs, even at temperatures T far above 
the bulk I-N transition temperature T^^v- '^^^ symmetry 
breaking doesn't occur spontaneously, as characteristic 
of a genuine phase transition, but is enforced over rele- 
vant distances by the interaction with the walls Q . Thus 
confinement here plays a similar role as an external mag- 
netic field for a spin system: The strong first order I-N 
transition is replaced by a weak first order or continuous 
paranematic-to-nematic {P-N) transition, depending on 
the strength of the surface orientational field. 

Whereas a qualitative understanding of this behavior 
has been achieved for a variety of spatially mesoconfined 
LCs P, more detailed comparisons with theoretical 
predictions have been extremely challenging in the past, 
mainly due to the complex, tortuous, multiply connected 



pore networks or unknown surface/interface-LC interac- 
tions in the aforementioned studies. The advent of arrays 
of straight nanochannels of silica and silicon with simpler 
channel geometries may allow to gain deeper, quantita- 
tive insights into this phenomenology. 

In this Letter, we present a high-resolution optical bire- 
fringence study of rod-hke nematogens (7CB and 8CB) 
confined to an array of parallel, non-tortuous channels 
of 10 nm mean diameter and 300 micrometer length in 
a monolithic silica membrane. We demonstrate that the 
optical transparency, the straight channel geometry of 
the host along with the sensitivity of the modulated beam 
technique employed allows us to precisely characterize 
the orientational order of the nematogens with respect 
to the long axis of the nanochannels. The simple type of 
restricted geometry will allow us to compare our measure- 
ments quantitatively with a Landau-de Gennes model for 
the P-N transition and, consequently, to determine both 
the strength of the orientational ordering field, imposed 
by the silica walls, and the influence of quenched disor- 
der, attributable to channel irregularities. 

For rod-like molecules the degree of orientational 
molecular order can be quantified by the uniaxial order 
parameter Q = ^ ^3cos^ 0—1), where (/) is the angle be- 
tween the long axis of a single molecule and a direction 
of preferred orientation of that axis, the director. The 
brackets denote an averaging over all molecules under 
consideration. The orientation of the director may vary 
locally. However, it can be dictated by external fields 
or by surface anchoring conditions over macroscopic dis- 
tances. Planar silica surfaces enforce planar anchoringof 
7CB and 8CB without a preferred lateral orientation [91]. 
Additionally, the director is expected to be oriented par- 
allel to the long axis in a cylindrical silica channel [lO|. 
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FIG. 1: (color online). Birefringence of 7CB and 8CB measured in the bulk state, panel (a), (b), and in the silica nanochannels, 
panel (c) and (d), resp., as a function of temperature in comparison to fits (solid lines) based on the KKLZ-model discussed 
in the text. The final birefringence characteristic of the paranematic phases are shaded down to the P — N "transition" 
temperatures, TpN. The dashed lines mark the bulk I-N and N-SmA transition temperatures. As insets in (a) and (c), the 
bulk isotropic (I^) as well as the bulk nematic (A^'') phases upon homeotropic alignment, and the confined paranematic (P) 
and nematic (N) phases are illustrated, respectively. 



A statement which shall be explored in the following by 
birefringence measurements. 

The propagation speed of light and thus the refractive 
index n in a LC sensitively depends on the orientation 
of the polarization with respect to the molecular orien- 
tation of the anisotropic nematogens. Conversely, the 
state of molecular order in a LC can be inferred from 
optical polarization measurements. To a good approx- 
imation, Q is proportional to the optical birefringence 
An = ric — rio, where Hq and Ue refer to polarizations 
perpendicular and parallel to the local optical axis, the 
so-called ordinary and extraordinary refractive indices 
ll|, respectively. In a nematic LC the local optical 
axis agrees with the director. Thus, in principle it is 
sufficient to determine the experimentally accessible An 
in order to determine the molecular arrangement in an 
LC. However, there are weak, but final T-dependencies 
of the bare refractive indices, which do not originate from 
changes of the averaged collective molecular orientations 
but from changes in the anisotropic molecular polariz- 
abilities of the single molecules as a function of T. In 
order to separate out these effects, we resort to the quan- 
tity An* = nl{T) - nl{T) oc An(T) • {n,{T) + no(T))/2 
which can be shown to be solely proportional to Q{T) 
fist . For simplicity, we will refer to An* as "effective 
birefringence" in the following. 

For our measurements, a monolithic silica mem- 
brane permeated by an array of parallel aligned, non- 



interconnected channels of 300 /im length was prepared 
by thermal oxidation of a free-standing silicon membrane 
hj at 800 °C for 12 hours. The mean channel diameter 
D and porosity P was determined by recording a vol- 
umetric nitrogen sorption isotherm at T = 17 K to 
Li^lO.OiO.S nm and P = 50±2%^espectively. Elec- 
tron micrographs of the channels [I4I indicate sizeable 
1.0±0.5 nm mean square deviations of their surfaces from 
an ideal cylindrical form. The membrane was completely 
filled with the LCs by spontaneous imbibition [l3|. For 
the bulk measurement, we used a 50 ^m thick glass cell 
containing homeotropically aligned LCs, see Fig. [Ija). 

We used a high-resolution optical modulated beam 
method for the accurate determination of the phase re- 
tardation R between two perpendicularly polarized com- 
ponents of light transmitted through the samples. The 
setup, see Fig. [51 employs an optical photoelastic modu- 
lator and a dual lock-in detection scheme in order to min- 
imize the influence of uncontrolled light-intensity fluctu- 
ations 0, [3| • After passing the sample the laser light in- 
tensity (A = 632.8 nm) was detected by a photodiode and 
two lock-in amplifiers, which simultaneously determined 
the amplitudes of the first (C/n) and second {U2Q) har- 
monics, respectively. The phase retardation by the sam- 
ple R = arctan[([/aJ2(^o))/(C/2nJi(Ao))] (here Ji(Ao) 
and ^2(^0) are the Bessel functions corresponding to the 
PEM retardation amplitude ^0 = 0.383A) was measured 
for an incident angle 6 = 43.5 deg between laser beam 
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and sample surface and thus between beam and long axes 
of the silica channels, see inset of Fig. [2l For such a 
tilted sample geometry the conversion of the retardation 
R to An and An* was performed by numerically solving 
Berek's compensator formula (lit . 

In Fig. [2 An* is plotted for bulk 7CB and 8CB 
upon slow cooling (~0.01 K/min) to the solidification 
temperature. There is a jump in An*{T) of bulk 
7CB typical of the first-order I-N phase transition at 
Tj'^N ~ 42 °C^. An*(r) of 8CB exhibits the signa- 
tures characteristic of the first-order I-N transition at 
^IN ~ 41 °C and an almost continuous N-SmA transi- 



tion at TNgniA °C[l|,ll^. Any pretransitional effects 
are clearly absent in the bulk isotropic phase of both LCs 
investigated, An*(r) ^ for T > T^^. 

The nanoconfined LCs reveal a considerably different 
behavior, see Fig. [T] lower panels. In agreement with 
neutron diffraction experiments on 8CB in nanochannels, 
there is no indication of a sharp N-SniA phase transi- 
tion [l3]. More interestingly, there exists a residual An* 
characteristic of a paranematic LC state at Ts far above 
Tj*^. Upon further cooling An* increases continuously 
and at the lowest Ts investigated, the absolute magni- 
tude of An* is compatible with an 80% (75%) alignment 
of the 7CBs' (8CBs') long axes parallel to the channel 
axes. Hence, the silica nanochannel confinement dictates 
indeed a substantial molecular alignment, as proposed in 
the introduction and illustrated in the inset of Fig. [TJ 
More importantly, it renders the transition continuous. 
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FIG. 2: (color online). Schematic experimental setup of a 
high-resolution birefringence measurement consisting of a He- 
Ne laser, an optical polarizer (P), a temperature controlled 
(TC) optical cell (OC), an optical analyzer (A), a photoelastic 
modulator (PEM-90), a photodiode as a detector (PD) and a 
"lock-in" detection and analyzing unit. The inset depicts the 
orientation of the silicon nanochannel membrane with respect 
to the incident laser beam. 

In the following we are going to analyze the peculiar 
An* behavior within a Landau-de Gennes model for the 
I — N transition in confinement suggested by Kutnjak, 
Kralj, Lahajnar, and Zumer (KKLZ-model) [2(, |6|]. The 
dimensionless free energy density of a nematic phase spa- 
tially confined in a cylindrical geometry with planar an- 
choring conditions reads in the KKLZ-model as: 



/ = ig2 _ + q'^ -qcr + nq^ 



(1) 



where q = Q/Q{Ti^) is the scaled nematic order param- 
eter, t is a reduced temperature, and a is the effective 
surface field. The last term in Eq. 1 describes quenched 
disordering effects due to surface-induced deformations 
(wall irregularities) [l[. Minimalization of / yields the 
equilibrium order parameter q^, which is shown in Fig. [3] 
for selected values of a and function of t. In 

the KKLZ-model, the I-N transition is of first order for 
a < (Jc = 0.5. The jump of qe approaches zero while 
a — > 0.5, see inset in Fig. [3l Thus, CTc marks a critical 
threshold separating first-order, discontinous from con- 
tinuous I-N behavior. 

In the following we apply the KKLZ-model to our 
measured An*(T). The solid lines in Figs. [IJa),(b) are 
the best fits of the dependencies An* (T) as obtained by 
rescaling qc and t while assuming an absence of any sur- 
face ordering and quenched disorder fields in the bulk 
state {(j{bulk) = a{D = C!o)=0, n = 0). Thereby, we 
achieve an encouraging agreement between the measured 
An* (T) curves in the proximity of T^^ and deep into the 
nematic phase for both bulk LCs. The sizeable deviations 
for 7CB below 230 °C originate in the neglect of higher 
order terms in the KKLZ-model, necessary to produce 
the saturation behavior of q at low T. 

More importantly, we achieve also an excellent agree- 
ment between measured and KKLZ- modelled An*(T) 
over the entire P-N transition regime for the nanocon- 
fined LCs, provided we assume final surface ordering 
fields of magnitude a{D = lOnm) = 0.81, 1.15 and 
strengths of the quenched disorder parameter k = 0.83, 
1.4 for 7CB and 8CB, resp., see panel (c) and (d) in 
Fig. [U respectively. As expected from the observed 
continuous behavior of An*(T) in both cases <j{D) is 
well above Cc = 0.5. Moreover, the upward shift in 
TpN due to the anchoring surface field, predicted by 
the KKLZ-model, is slightly overbalanced by a down- 
ward shift due to quenched disorder effects for this set of 
KKLZ-parameters - see also Fig. |31 

The bare geometrical mechanism which favors an align- 
ment of rod- like molecules in cylindrical confining geome- 
tries is expected to increase in strength with the length 
of the molecules [l^l- Also the alignment effect of the 
anchoring field of the confining walls increases with this 
length. Therefore the cr-increase of 30% appears not 
too surprising, despite the relatively small 10% length- 
increase between 8CB and 7CB. 

According to the KKLZ-model cr{D) scales with 1/D, 
which also allows us to estimate cr as a function of D for 
other confining silica geometries. In particular, we arrive 
at a Dc of 23 nm and 16 nm for the "critical" silica pore 
diameter of 8CB and 7CB, resp., {a{Dc) — Uc — 0.5) sep- 
arating continuous from discontinuous behavior. This is 
in quantitative agreement with the still weakly discontin- 
uous behavior reported for 8CB confined in 24 nm xerogel 
mesopores f2| and the continuous behavior found for 7CB 
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in 7 nm mean pore diameter Vycor glass [l| 




FIG. 3: (color online). Order parameter qc in the KKLZ- 
model as a function of reduced temperature t for selected 
effective surface fields a in the absence of quenched disorder, 
n — (solid lines) and for the a-, k- values which yield the best 
fits of An*{T) measured for the confined LCs (dashed lines). 
Inset: Order parameter jump at the I-N (P-N) transition as 
a function of the surface field a (k = 0). 

In conclusion, our measurements indicate that in 
10 nm straight silica channels the surface anchoring fields 
render the bulk discontinuous I-N transition to a con- 
tinuous P-N transition. The transition temperature is 
changed only marginally, due to a balance of its surface 
ordering induced upward and its quenched disorder in- 
duced downward shift, similarly as has been observed 
for LCs imbibed in tortuous Vycor glass pores and a 
variety of other mesoporous matrices The purely 
phenomenological findings presented here would profit 
from more microscopic information, for example gained 
by Monte Carlo or Molecular Dynamics simulations [7,] 
or from x-ray and neutron scattering experiments. 

Finally we would like to state, that the rheology of 
bulk LC changes abruptly at the I-N transition [19, 20]. 
By contrast, the gradual P-N transition reported here 
should lead to a continuous T-evolution of the fluidity of 
LC in silica nanochannels. Moreover, velocity slippage 
at the walls, crucial in the emerging field of nanofiuidics 
is expected to be associated with the molecular 
alignments in the channels [l^l- Thus, we hope the pecu- 
liarities reported here will stimulate LC flow experiments 
with silica nanochannels, which would be important ex- 
tensions of previous ex per iments on the fluidity of LCs 
in thin film geometries [23| . 
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